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Abstract 

We calculate three-loop corrections to correlation functions of heavy-quark currents 
in the low- and high-energy regions. We present 30 coefficients both in the low-energy 
and the high-energy expansion of the scalar and the vector correlator with non- 
diagonal flavour structure. In addition we compute 30 coefficients in the high-energy 
expansion of the diagonal vector, axial-vector, scalar and pseudo-scalar correlators. 
Possible applications of our new results are improvements of lattice-based quark- 
mass determinations and the approximate reconstruction of the full momentum 
dependence of the correlators. 

Key words: Perturbative calculations, Quantum Chromodynamics, Heavy Quarks 
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1 Introduction 

The precise determination of the fundamental parameters of the Standard 
Model (SM) is of utmost importance to pin down the fundamental interac- 
tions of the elementary particles and to establish deviations from the SM. 
Three especially interesting parameters, namely the strong coupling constant 
and the masses of charm and bottom quarks, can be determined with high 
accuracy from correlators of heavy-quark currents by comparing their mo- 
ments with experimental data using sum rules jTP] . Conventionally, moments 
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from the low-energy expansion have been used for this purpose j3llHl5|6f7f8] . 
but also methods relying on high-energy moments have been shown to yield 
quite accurate results [9j. 



As an alternative to experimental data also lattice simulations can be used 
in order to compute moments of two-point functions. Thanks to significant 
improvements in the last years results can be obtained not only for the charm- 
quark mass and the strong coupling constant [IQl, but also for the bottom- 
quark mass [llj. 

Aside from these determinations, quark-current correlators have a number of 
additional phenomenological applications. Two examples are the sum-rule de- 
terminations of heavy-meson decay constants and as{mr) (for recent works 
see e.g. Refs p^l3pi4|[T5] ). While the former requires - among other ingre- 
dients - the low-energy expansion of the respective correlator, the latter uses 
the absorptive part above the production threshold. 

In this work we consider correlators of quark currents with various different 
Lorentz structures. Furthermore we distinguish between diagonal and non- 
diagonal correlators. The diagonal correlators are defined through 

n§(?) = ^/e^''^(0|J^(x)J^(0)|0)da;, (1) 

with the currents J^{x) given by 

J^{x) =ip{x)T^^{x). (2) 
The non-diagonal correlators are defined in analogy by 

^Nuiq) = ^ I e^^-(0|j^(a;)j^(0)|0)da:, (3) 
with the currents j^{x) given by 

fix) = V'(x)r^x(a:) . (4) 

Here ip denotes a heavy quark while x stands for a massless one. The scalar, 
pseudo-scalar, vector and axial-vector currents are defined through the Dirac 
structures T-^ = 1, ^75, 7^^, 7/^75 {X = s,p,v,a), respectively. Note that in the 
case of the non-diagonal correlators Hf^D = IInd ^^^1 n^D = hold. 

For the determination of the quark masses from experimental data only the 
diagonal vector correlator 11^ can be used while it is in principle possible 
to use all diagonal and non-diagonal correlators in combination with lattice 
simulations. 

Current-current correlators have already been extensively studied in the li- 
terature. The two-loop corrections to the diagonal vector correlator have been 
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calculated in [16], the corresponding corrections for the non-diagonal one 
in [T7j . At three loops the leading terms of the low- and high-energy expansion 
have been determined in jTHP^ . For the diagonal correlators many terms in 
the low-energy expansion have been computed in [20II2T] . At four loops the 
leading terms in the low-energy expansion of various currents have been cal- 
culated in [22][23]IHl2m25] ■ In the high-energy region the expansion at four- loop 
order has been obtained in [26]. In [2^ effects of a second massive quark have 
been taken into account at three-loop order. Collecting the available informa- 
tion for the correlators at three and four loops, the full momentum dependence 
has been approximated using Fade approximations and Mellin-Barnes inspired 
methods [28f 18f29ll30f31f32j . In this paper we complete the calculation of the 
low- and high-energy expansion of diagonal and non-diagonal heavy-quark 
correlators using an improved method for the calculation. 

In Section |2] we present the details of the calculation. In Section |3] we present 
explicit results for the non-diagonal vector correlator in numerical form and 
discuss the convergence of the expansions. The analytical results for all corre- 
lators can be found on the web page 

http://www-ttp. particle. uni-karlsruhe.de/Frogdata/ttpll/ttpll-25/[ 



2 Details of the calculation 



The calculation is organized as follows. The diagrams are generated using 
qgraf [33] and are subsequently mapped onto a small set of 11 topologies 
(cf. Fig. [T]) using q2e and exp [M]. After the application of suitable projec- 
tors the resulting scalar integrals are reduced to the master integrals shown 
in Fig. [2] using integration-by-parts identities [351136] implemented in the pro- 
gram Crusher [37J. The resulting expression which still contains the full de- 
pendence on both the external momentum square and the heavy-quark 
mass m is then - lacking an analytical result for the master integrals - ex- 
panded in the energy region of interest as explained in more detail below. 
The expanded master integrals are combined in a FORM [38] program in order 
to obtain the final expressions valid in either the low-energy or the high- 
energy region. For further use the low- and high-energy expansions of all mas- 



ter integrals can be obtained from the web site http://www-ttp.particle.uni 



karlsruhe .de/Progdata/ttpll/ttpll-25/ masters / 



The coupling constant is renormalized in the MS scheme while for the 
renormalization of the heavy-quark mass m we employ both the MS and the 
on-shell scheme. In the following we concentrate on the renormalization in the 
MS scheme, the results for the on-shell scheme can be found on the web page 
mentioned above. For the overall normalization we employ 11"^ (0) = with 
X G a, s,p. 
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Fig. 1. Three- loop master topologies needed for the calculation of the non-diagonal 
and diagonal correlators. Solid lines denote massive propagators while dotted lines 
represent massless ones. 

Let us now have a closer look at the technical details of the expansion of the 
master integrals. The starting point for the calculation of the expansion is a 
suitable differential equation, which can easily be obtained from the master 
formula [SMUT] 



' 2 5 2 ^ ^ 
dq^ dm? 



(5) 



The derivative 



where is the mass dimension of the integral Mi{q^,m? 
with respect to m? can be performed resulting in a differential equation with 
respect to q^ for the integral Mi{q'^, m?). The expansion in 2; = q^ /m? can now 
be obtained by inserting a suitable ansatz, e.g. 



2 2^, 

,m J 



k=0 



(6) 



into the differential equation and solving the resulting linear system of equa- 
tions for the coefficients fc. Ansatz ^ is sufficient to obtain the low-energy 
expansion for the non-diagonal correlators. The only further input needed is 
the value of the integral for z = 0, which corresponds to the initial condition 
of the differential equation. In the case of the high-energy expansion ansatz 
(|6| is not sufficienip"! but one has to use an ansatz of the form 



2e) 
k I 



k=k{\ 



instead. The lower limit and the initial conditions Aj^^g, A-'^^, / 
can be determined from the leading terms of the asymptotic expansion for 
z — )■ —00. Generally, these leading terms consist of products of massive vacuum 
diagrams and massless propagators with a total number of three loops in 

^ An extension is also necessary for the low-energy expansion of the diagonal cor- 
relators as discussed in detail in Ref. I21i. 



(7) 
1,2,3 
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Fig. 2. All master integrals needed for the calculation of the diagonal and non-di- 
agonal correlators up to three loops. Solid lines denote massive propagators while 
dotted lines represent massless ones. A dot on a line denotes a squared propagator. 

each product. The coefficients ^f^^ are given by the products of propagators 
with a total of I loops and vacuum diagrams with 3 — / loops. Consequently, 
Ai^ko consists of the products without any massless propagators. The diagrams 
contributing to the initial conditions are shown in Fig. [3] in order to obtain 
dimensionless coefficients Aj ^g, Af^^^ we set — = = 1 in those diagrams. 

It should be noted that in practice it is convenient to just use a small enough 
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universal value for ko instead of carefully deriving it separately for each single 
master integral from the respective asymptotic expansion. Furthermore, most 
of the coefficients Ai^^o^ ^f^ko already determined by the linear system 
of equations and only a few remaining ones have to be determined through 
asymptotic expansion. The considerable freedom in the choice of the initial 
conditions also implies that the set of contributing diagrams shown in Fig. [3] 
is by no means unique. 

Finally, let us show as an example the explicit high-energy expansion of one 
specific master integral: 

M6l(g^m2) 

+ {3z-' + z-^ + ... . 

-48z-2 + 5;z-^ + ...)--- 

+ f4z-=^ + 2;z-2 + ...V-<' '■>--'( 

Af\-^ =( - mt-\-^ + 36e-2^-2 - Qt-^z-^ + bUt'^z-^ - 2Mt-^z-^ 
+ 29e~^;z"^ - 2665^"^ + 508^"^ - AQz'^ 
+ (-3^-3-z-^-^-i + ...)-<(T>--. (8) 

In this example terms of order e or 1/ z'^ and higher have been omitted in the 
coefficients of the integrals. 

In the case of the pseudo-scalar and axial-vector correlators one has to ad- 
dress the treatment of 75. In this calculation we use the prescription given by 
Larin [l2] for the so-called singlet diagrams with exactly one 75 in a fermion 
trace and a naively anticommuting 75 for all other diagrams. In order to 
present only anomaly-free quantities the axial-vector correlator is computed 
for a quark isospin doublet {ip, x) like in Ref. j21| . 





3 Results 



Since the full results are very lengthy we present only the numerical results 
for the non-diagonal vector correlator. The analytical results for all correlators 
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Fig. 3. Building blocks for the initial conditions ^i,fcoi ^ili' ^ ~ 1)2,3. The initial 
conditions are linear combinations of products of the diagrams shown here. The 
diagrams are evaluated for — = m? = 1 in order to arrive at dimensionless 
coefficients Ai^^^, Afj}^. 

can be found on the web page 

[http://www-ttp.particle.uni-karlsruhe.de/Progdata/ttpll/ttpll-25/ 

The expansion has been performed in 2; = /m?. The results for the first 6 
low-energy and 8 high-energy coefficients agree with Ref. [T9] . 

The vector correlator can be written as 

n^D,.!?') = (-^'^7"^ + q'q'WUq') + , (9) 

We present only the results for the transversal part of the vector correlator, 
the longitudinal part can be obtained from the results for the scalar one. In 
the following we set the renormalization scale n = m during the numerical 
evaluation. 

In Table [T] we present the numerical results for the low-energy expansion of 
the non-diagonal vector correlator which is decomposed according to 




C^^''"[nh\ and C^^''"[ni] contain the contributions from closed heavy and light 
quark loops, respectively. 

The results for the high-energy expansion shown in Tables 2-4 are more com- 
plicated due to logarithmic contributions. Since the correlator is known ana- 
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lytically up to two loops we present only the results for the three-loop contri- 
bution. It is decomposed as 



n 



i), 3— loop 
ND 



n>0 m=0 



f)m,v 

n,m 



ni z 



(11) 



In Fig. |4]we show the convergence of the series for the case of the non-diagonal 
vector correlator. The behavior of the series is shown both for the low-energy 
and the high-energy expansion. For the high-energy expansion we show the 
real and the imaginary part. Since the non-diagonal correlators have a cut 
through three heavy quark lines, the high-energy expansion converges only 
above z = 9, which corresponds to a velocity of the heavy quark of f = 
{z — l)/{z + 1) = O.Sj^This behavior can also clearly be observed in Tables 
2 and 3 which show rapidly growing coefficients corresponding to a radius of 
convergence r < 1. Since Table 4 contains only contributions from diagrams 
with a closed massless- quark loop, the coefficients are better behaved. 



^ The diagonal correlators have a four-particle cut starting at three-loops and con- 
verge therefore only above z = 16, i.e. v = y^l — 1/z = 0.97. 
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n Cf'^ (7(^)'^ C^^^'ll] CF'^'K] C^^^'''[ni] 



1 


0, 


.500000 


0.457734 




1.51195 


0.0166564 


0, 


.337837 


2 


0, 


.133333 


-0, 


.0483080 




0.840974 


-0, 


.0369934 


0, 


.187255 


3 


0, 


.0555556 


-0, 


.115602 




0.665554 


-0, 


.0289884 


0, 


.135615 


4 


0, 


.0285714 


-0, 


.115219 




0.475376 


-0, 


.0213752 


0, 


.104474 


5 


0, 


.0166667 


-0, 


.102152 




0.309702 


-0, 


.0161172 


0, 


.0832408 


6 


0, 


.0105820 


-0, 


.0880811 




0.177592 


-0, 


.0125094 


0, 


.0679630 


7 


0, 


.00714286 


-0, 


.0756354 




0.0753998 


-0, 


.00996452 


0, 


.0565728 


8 


0, 


.00505051 


-0, 


.0651894 


0, 


.00291844 


-0, 


.00811379 


0, 


.0478459 


9 


0, 


.00370370 


-0, 


.0565463 


0, 


.0628596 


-0, 


.00672995 


0, 


.0410083 


10 


0, 


.00279720 


-0, 


.0494003 


0, 


.108797 


-0, 


.00566985 


0, 


.0355496 


11 


0, 


.00216450 


-0, 


.0434648 


0, 


.144058 


-0, 


.00484063 


0, 


.0311210 


12 


0, 


.00170940 


-0, 


.0385013 


0, 


.171144 


-0, 


.00418016 


0, 


.0274778 


13 


0, 


.00137363 


-0, 


.0343196 


0, 


.191926 


-0, 


.00364579 


0, 


.0244440 


14 


0, 


.00112045 


-0, 


.0307699 


0, 


.207816 


-0, 


.00320745 


0, 


.0218902 


15 


0, 


.000925926 


-0, 


.0277347 


0, 


.219882 


-0, 


.00284350 


0, 


.0197199 


16 


0, 


.000773994 


-0, 


.0251214 


0, 


.228940 


-0, 


.00253805 


0, 


.0178596 


17 


0, 


.000653595 


-0, 


.0228568 


0, 


.235620 


-0, 


.00227923 


0, 


.0162527 


18 


0, 


.000556948 


-0, 


.0208824 


0, 


.240408 


-0, 


.00205802 


0, 


.0148551 


19 


0, 


.000478469 


-0, 


.0191512 


0, 


.243686 


-0, 


.00186749 


0, 


.0136316 


20 


0, 


.000414079 


-0, 


.0176253 


0, 


.245754 


-0, 


.00170221 


0, 


.0125545 


21 


0, 


.000360750 


-0, 


.0162738 


0, 


.246850 


-0, 


.00155793 


0, 


.0116011 


22 


0, 


.000316206 


-0, 


.0150713 


0, 


.247164 


-0, 


.00143123 


0, 


.0107532 


23 


0, 


.000278707 


-0, 


.0139968 


0, 


.246847 


-0, 


.00131937 


0, 


.00999562 


24 


0, 


.000246914 


-0, 


.0130328 


0, 


.246023 


-0, 


.00122012 


0, 


.00931596 


25 


0, 


.000219780 


-0. 


.0121648 


0, 


.244790 


-0, 


.00113166 


0, 


.00870382 


26 


0, 


.000196483 


-0, 


.0113806 


0, 


.243229 


-0, 


.00105247 


0, 


.00815053 


27 


0, 


.000176367 


-0, 


.0106696 


0, 


.241406 


-0, 


.000981314 


0, 


.00764873 


28 


0, 


.000158907 


-0, 


.0100231 


0, 


.239376 


-0, 


.000917133 


0, 


.00719220 


29 


0, 


.000143678 


-0, 


.00943357 


0, 


.237182 


-0, 


.000859046 


0, 


.00677564 


30 


0, 


.000130336 


-0, 


.00889448 


0, 


.234861 


-0, 


.000806306 


0, 


.00639448 



Table 1 

Numerical results for the low-energy moments for the non-diagonal vector correlator 
in the MS scheme. 
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n 


71, U L J 


71,1 L J 


71, Z L J 


n,o L J 





1.78284 


-2.64761 


1.83333 





1 


-17.7072 


25.1337 


-16.3333 


3.16667 


2 


-60.1359 


47.0000 


-8.00000 





3 


-41.3176 


-15.4431 


16.0772 


-3.32442 


4 


4.27239 


-20.8764 


5.72599 


0.659122 


5 


7.18683 


-11.8251 


0.822323 


-0.129739 


6 


0.493132 


-5.22086 


2.92024 


-0.566886 


7 


-16.8712 


-10.5240 


9.71602 


-0.991008 


8 


-56.8660 


-43.6415 


29.4833 


-1.70280 


9 


-164.387 


-190.388 


94.9042 


-2.87364 


10 


-485.549 


-831.419 


339.294 


-4.95108 


11 


-1537.17 


-3780.76 


1357.43 


-8.63056 


12 


-5240.01 


-18141.8 


6001.46 


-15.2756 


13 


-18953.8 


-91894.0 


28790.4 


-27.3450 


14 


-70903.6 


-488762 


147374 


-49.4690 


15 


-264436 


-2.71202 • 10^ 


794591 


-90.2629 


16 


-922552 


-1.56052 • 10'^ 


4.46970 • 10^ 


-165.935 


17 


-2.52510 • 10^ 


-9.26482 • 10'^ 


2.60485 • 10^ 


-306.988 


18 


-142055 


-5.65199 • 10^ 


1.56448 • 10^ 


-571.093 


19 


8.11996 • 10^ 


-3.53100 • 10^ 


9.64408 ■ 10^ 


-1067.52 


20 


9.74537 • 10^ 


-2.25279 • 10^° 


6.08196 ■ 10^ 


-2003.91 


21 


8.94750 • 10^ 


-1.46443 ■ 10^1 


3.91351 • 10^° 


-3775.71 


22 


7.44796 ■ 10^0 


-9.68039 ■ 10^1 


2.56374 ■ 10^1 


-7137.78 


23 


5.91860 • 10^^ 


-6.49650 • 10^2 


1.70673 • 10^2 


-13533.8 


24 


4.59022 • 10^2 


-4.41990 • 10^3 


1.15282 • 10^3 


-25729.9 


9^ 
ZO 


O.0i«j4'i • iU 


q n/i/isn 1 ni4 




— 4iyuo0.o 


26 


2.67074 • 10^'' 


-2.12156 • 10^^ 


5.46505 • lO^"^ 


-93657.0 


27 


2.02372 • 10^5 


-1.49382 • 10^6 


3.82727 • 10^^ 


-179244 


28 


1.53213 • 10^^ 


-1.06201 • 10^^ 


2.70750 • lO^*^ 


-343676 


29 


1.16067- 10^^ 


-7.61776 • W^'^ 


1.93326 • 10^^ 


-660072 


30 


8.80664 • 10^'^ 


-5.50954 • 10^^ 


1.39237 • 10^^ 


-1.26974-10 



Table 2 

Three-loop coefficients of the high-energy expansion for the non-diagonal vector 
correlator in the MS scheme with = nj = 0. 
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n 




71,1 L '(-J 









-0.251201 


0.153727 


-0.111111 





1 


0.154182 


-1.14726 


0.555556 


-0.111111 


2 


3.61746 


-1.14459 


-0.666667 





3 


1.27971 


2.67597 


1.06790 


0.185185 


4 


-3.99022 


-1.81636 


-0.861111 


-0.543210 


5 


-8.59452 


0.842895 


1.23685 


-1.15062 


6 


-25.2903 


-2.26222 


8.70815 


-1.68560 


7 


-65.3844 


-31.4942 


28.2131 


-2.18624 


8 


-164.223 


-160.905 


86.7525 


-2.66878 


9 


-441.129 


-715.641 


295.743 


-3.14070 


10 


-1320.45 


-3233.48 


1154.83 


-3.60588 


11 


-4384.76 


-15440.6 


5072.08 


-4.06652 


12 


-15690.9 


-78173.7 


24362.3 


-4.52396 


13 


-58311.1 


-416730 


125176 


-4.97906 


14 


-215239 


-2.32045 • 10^ 


677824 


-5.43239 


15 


-732644 


-1.34050-10^ 


3.82947 - 10^ 


-5.88438 


16 


-1.83932 • 10^ 


-7.99084 - 10^ 


2.24126 - 10'^ 


-6.33529 


17 


1.84847 • 10^ 


-4.89431 - 10^ 


1.35168 - 10^ 


-6.78534 


18 


8.23216 • 10'^ 


-3.06952 - 10^ 


8.36569 - 10« 


-7.23470 


19 


9.21471 • 10^ 


-1.96567- 10^° 


5.29604 - 10^ 


-7.68348 


20 


8.27344 • lO'^ 


-1.28233-10^^ 


3.42033 - 10^° 


-8.13178 


21 


6.81847 ■ 10^° 


-8.50550 - 10^1 


2.24852 - 10^1 


-8.57968 


22 


5.39132 • 10^1 


-5.72647 - 10^2 


1.50188 - 10^2 


-9.02723 


23 


4.17068-10^2 


-3.90796 - 10^3 


1.01766-10^3 


-9.47449 


24 


3.18860-10^3 


-2.69998 - 10^^ 


6.98579 - 10^3 


-9.92149 


25 


2.42300 - 10^^ 


-1.88650 - lO^^ 


4.85260 - 10^^ 


-10.3683 


26 


1.83633 - lO^s 


-1.33179 - 10^6 


3.40755 - 10^5 


-10.8149 


27 


1.39098 - 10^6 


-9.49177-10^6 


2.41677- 10^6 


-11.2613 


28 


1.05453 - lO^'^ 


-6.82454 - lO^'^ 


1.72987- 10^'^ 


-11.7075 


29 


8.00841 - 10^^ 


-4.94695 - 10^^ 


1.24877- lO^s 


-12.1537 


30 


6.09583 - 10^8 


-3.61318-10^^ 


9.08598 - lO^s 


-12.5997 



Table 3 

Three-loop contribution proportional to Uh to the coefficients of the high-energy 
expansion for the non-diagonal vector correlator in the MS scheme. 
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n 


71, U L 


71, i L 'J 




7Z,o L tJ 





-0.730175 


0.153727 


-0.111111 





1 


-0.314922 


-1.14726 


0.555556 


-0.111111 


2 


0.984271 


-1.77778 








3 


1.27530 


-0.0592818 


-0.615226 


0.0534979 


4 


0.544450 


0.835905 


-0.0709877 


-0.0246914 


5 


0.0510031 


0.410862 


0.0306790 


-0.00740741 


6 


-0.0437552 


0.209240 


0.0303704 


-0.00329218 


7 


-0.0561161 


0.117605 


0.0233371 


-0.00176367 


8 


-0.0503344 


0.0722787 


0.0177016 


-0.00105820 


9 


-0.0420848 


0.0475365 


0.0136697 


-0.000685871 


10 


-0.0347091 


0.0329283 


0.0107911 


-0.000470312 


11 


-0.0287111 


0.0237565 


0.00869670 


-0.000336700 


12 


-0.0239493 


0.0177088 


0.00713792 


-0.000249408 


13 


-0.0201760 


0.0135592 


0.00595223 


-0.000189934 


14 


-0.0171671 


0.0106166 


0.00503229 


-0.000148000 


15 


-0.0147456 


0.00847129 


0.00430581 


-0.000117578 


16 


-0.0127769 


0.00687002 


0.00372304 


-0.0000949668 


17 


-0.0111601 


0.00565028 


0.00324902 


-0.0000778089 


18 


-0.00981938 


0.00470458 


0.00285863 


-0.0000645526 


19 


-0.00869741 


0.00395987 


0.00253355 


-0.0000541477 


20 


-0.00775052 


0.00336532 


0.00226014 


-0.0000458663 


21 


-0.00694509 


0.00288479 


0.00202813 


-0.0000391926 


22 


-0.00625501 


0.00249212 


0.00182965 


-0.0000337544 


23 


-0.00565976 


0.00216805 


0.00165859 


-0.0000292783 


24 


-0.00514313 


0.00189817 


0.00151016 


-0.0000255604 


25 


-0.00469214 


0.00167159 


0.00138059 


-0.0000224467 


26 


-0.00429635 


0.00147992 


0.00126683 


-0.0000198192 


27 


-0.00394727 


0.00131666 


0.00116643 


-0.0000175864 


28 


-0.00363797 


0.00117672 


0.00107740 


-0.0000156771 


29 


-0.00336274 


0.00105607 


0.000998082 


-0.0000140345 


30 


-0.00311683 


0.000951473 


0.000927135 


-0.0000126137 



Table 4 

Three-loop contribution proportional to ni to the coefficients of the high-energy 
expansion for the non-diagonal vector correlator in the MS scheme. 
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4 Conclusion 



We have used modern techniques to compute 30 terms in the low-energy ex- 
pansion of the non-diagonal scalar and vector correlators, matching the al- 
ready available information on diagonal correlators [21]. Using slightly modi- 
fied methods we have also obtained 30 coefficients in the high-energy expansion 
of both non-diagonal and diagonal correlators. Due to multiply massive cuts 
in some of the contributing diagrams it is expected that the high-energy ex- 
pansion already breaks down significantly above the physical threshold. Using 
only the previously available coefficients this behaviour is rather difficult to 
see. The newly computed coefficients, however, clearly support the presence 
of such a divergence. 

As a final remark it should be noted that the expansion via differential equa- 
tions used in this publication is rather efficient. This means that it is easy to 
obtain significantly more terms in the expansions if the need arises. 



Acknowledgments 

This work was supported by the Deutsche Forschungsgemeinschaft through 
the SFB/TR-9 "Computational Particle Physics". A.M. thanks the Landes- 
graduiertenforderung for support. We would like to thank T. Kasprzik, J.H. Kiihn 
and M. Steinhauser for carefully reading the manuscript. 



13 




-10 



N = Q 

N = 12 

AT = 18 

N = 2i 

AT = 30 



J I I I I I I I I I Llj L 



_L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

V 



10 

SI 
G 

S 



-| I I I I I I I I I l| I I I [T 



N = 6 
N = 12 
■N =18 
■N = 2A 
iV = 30 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

V 

Fig. 4. Convergence of the three-loop contribution to the low- and high-energy 
expansion of the non-diagonal vector correlator in the MS scheme for n/i = 1, 
ni = 3. The upper figure shows the correlator in the low-energy range for different 
orders of the expansion in z. The lower figures show the behavior above threshold 
for both real and imaginary parts. Note that the high-energy expansion converges 
only above v = 0.8 which corresponds to the three-particle cut. 
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